Surface impedance of superconductors with weak magnetic impurities 
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Magnetic impurities affect the spectrum of excitations of a superconductor and thus influence 
its impedance. We concentrate on the dissipative part of the surface impedance. We investigate 
its dependence on frequency, the density and strength of magnetic impurities, and the density and 
temperature of quasiparticles. Even a small concentration of weak magnetic impurities significantly 
modifies the excitation spectrum in the vicinity of the BCS gap. Therefore, we give special attention 
to the absorption threshold behavior at zero temperature and to the low-frequency absorption by 
quasiparticles. The discrete energy states introduced at low density of magnetic impurities may serve 
as traps for nonequilibrium quasiparticles, reducing the absorption in some range of low radiation 
frequencies. 

PACS numbers; 74.25.nn, 74.25.N-, 75.20.Hr, 73.20.Hb 
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I. INTRODUCTION 

Electron scattering off magnetic impurities, unlike the 
potential scattering, substantially modifies the properties 
of s-wave superconductors. It waSr^ealized in the semi- 
nal work by Abrikosov and Gor'kovid (AG) that the pres- 
ence of magnetic impurities reduces the superconducting 
gap below its BCS value and eventually may lead to the 
phenomenon of gapless superconductivity. The gap sup- 
pression was investigated by means of tuuDfiling between 
normal and superconducting electrodes.alj Beyond the 
AG theory, it was realized that a single magnetic wwu- 
rity creates a localized state within the BCS gap.EH3 A 
small concentration of magnetkrimpurities results in the 
formation of an impurity bandElfl which merges with the 
continuum as the concentration exceeds a certain value 
(see Ref. [sjfor a review). (The original AG theory is valid 
above that value.) The impurity band was investigated, 
by tunneling experiments with alloys (suck, as PbMn)El 
and normal metal-superconductor bilayers,E3 and also bj; 
means of thermal transport in superconducting films.til 
The observation of discrete levels associated with a sep- 
arate magnetic impurity came l ater , in the scanning tun- 
neling microscopy experiments. 

The effect of magnetic impurities came under scrutiny 
recently in connection with limits onj-jperformance 
of supercondueting elements in qubitsLa and pho- 
ton detectors.ta In these applications, superconduc- 
tors should perform in the ac regime. There is also 
some experimental effort directly aimed at observa- 
tion of the effect of magnetic impurities ©a the surface 
impedance of superconducting-multilayeralj and uncon- 
ventional supescoaductors.EIro Theoretically, the results 
on impedancell3Tlj were analyzed in the framework of a 
simple two-fluid model. The two-fluid model assumes a 
simplified frequency-independent Drude-type dissipative 
conductivity a, while the infiuence of magnetic impuri- 
ties on the electron spectrum should lead to a nontrivial 
a{u!) dependence on the frequency uj. 



A theoretical approach to studying the surface 
impedance in superconductors was developed in the cele- 
brated work of Mattis and Bardeen,E3 which treated the 
anomalous skin effect (the skin depth is smaller than the 
mean free path) taking nonmagnetic disorder into ac- 
count. Independently, the clean case .was, considered by 
Abrikosov, Gor'kov, and KhalatnikovOEj in the Green- 
function language. Later, NamE^ employed the same 
technique and demonstrated that the extreme anoma- 
lous limit is equivalent (up to the expressions for the 
normal-metallic conductivities) to the local limit, real- 
ized in the dirty regime. The complex conductivity of a 
superconductor containing both potential and magnetic 
impurities was addressed by Skalski et alE3 Their pa- 
per does contain general results applicable in the case 
of weak magnetic scattering; however, similarly to the 
AG theory, scattering off magnetic impurities was con- 
sidered in the Born approximation, which did not allow 
to treat the effect of localized states on the impedance. 
In Ref. |2^, Larkin and Ovchinnikov considered supercon- 
ductors with inhomogeneous pairing, which under cer- 
tain assumptions is formally equivalent to the presence 
of magnetic impurities in the AG regime. Some of their 
results can therefore be applied to this problem. How- 
ever, a convenient formulation of the theory, suitable for 
the experimentally relevant dirty limit and making it pos- 
sible to treat both the AG regime and the impurity band, 
as well as detailed quantitative investigation of the effect 
of magnetic impurities on the surface impedance in dif- 
ferent regimes over temperature, frequency, and electron 
scattering rate off magnetic impurities is still lacking, to 
the best of our knowledge. 

In this paper, we study the infiuence of magnetic impu- 
rities on the complex conductivity and surface impedance 
of a superconductor. We are especially interested in the 
limit of low frequency of the impinging radiation and low 
temperature, and in the effect of the impurity band on 
the dissipation. Quantitative results are presented for 
the experimentally relevant dirty limit (with respect to 
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the nonmagnetic disorder). 

The paper is organized as follows. First, in Sec. we 
consider the limit of weak scattering off magnetic impuri- 
ties and disregard nonmagnetic disorder ("clean" limit). 
In that limit, we find the localized states bound to mag- 
netic impurities and discuss the formation of the impu- 
rity band at finite but small concentration of magnetic 
impurities and the formation of the AG density of states 
at higher concentrations. We also present there a quali- 
tative consideration of the dissipative conductivity. The 
main features of the developed simple picture carry over 
to the realistic case of a "dirty" superconductor doped 
with magnetic impurities. In Sec. Ill, we consider the 
dirty limit for the nonmagnetic disorder in the framework 
of the quasiclassical formalism of the Usadel equation. 
We first obtain the spectrum of quasiparticle excitations 
in the presence of magnetic impurities (treated beyond 
the Born approximation). Then we address the effect of 
magnetic impurities on the low-temperature dissipative 
conductivity. In addition to the low- frequency limit, we 
also consider in detail the frequency dependence of the 
dissipation near the thresholds determined by the gaps 
in the excitation spectrum. Our results are summarized 
in Sec. IV. Some technical details are presented in the 
Appendixes. 



{ip and a are the operators of electrons and Bogoliubov 
quasiparticles, respectively): 

JS ^ crV-ki^-iAkaa 

ki,k2,o- 

= JS ^ [(wkiMk2 + WkiWk2)(alcit"k2T ~ "ki;"k2;) 



kl,k2 



-(wkiWk2 - "k2Wki)(aki^a^k2; ~ a-kiiak2r)J . (3) 

where ki(2) are momenta and a is the spin index (-1- 
for t, — for i). At the same time, the standard BCS 
part of the Hamiltonian becomes diagonal after the Bo- 
goliubov transformation, with the excitation energies 
i?k = \/$k + the superconducting state. Here A 

is the BCS gap and are the excitation energies in the 
normal state. Close to the gap (|i?k — A| ^ A), we can 
approximate E'k ~ A + ^k/2A, while the Bogoliubov co- 
herence factors are simply Uk ~ I'k ~ 1/ V^, so finally 
the total Hamiltonian takes the form 

H = ^ ekaL^ko- + JS ^ o-al^^ak^a, (4) 

k,CT ki,k2,cr 

where Ek = ^'^/2A and we measure the quasiparticle en- 
ergies from the edge of the continuous spectrum A: 



II. CLEAN LIMIT 



e = E-A. 



(5) 



To highlight the main features of the problem related 
to the impurity states, we consider here a simplified case 
of a superconductor with magnetic impurities only. In 
addition, considering the hybridization of single-impurity 
states into an impurity band we first assume magnetic 
impurities to be polarized (magnetic moments pointing 
up the z axis), while the impurity positions are random. 
Furthermore, we assume the exchange constant J of in- 
teraction of localized classical spins S with the spins of 
itinerant electrons to be small. Introducing the dimen- 
sionless parameter 



(1) 



describing the "strength" of a single magnetic impurity, 
we can write the latter requirement as 



C«i. 



(2) 



Here i/q = kp/2TT'^VF is the electron density of states 
(DoS) per spin projection at the Fermi level in the nor- 
mal state, with kp and vp being the Fermi wave vec- 
tor and the Fermi velocity, respectively. Weak exchange 
interaction will modify substantially the excitation spec- 
trum only close to the gap. Concentrating on that energy 
region, we simplify the Hamiltonian of the system in the 
following way. First, we express the operator of exchange 
interaction of electrons with an impurity (placed at ori- 
gin, R = 0) in terms of the Bogoliubov quasiparticles 



A. Single-impurity bound state 

It is easy to see that Hamiltonian (||) has a bound 
state for quasiparticles with spin a satisfying the condi- 
tion JSa < 0. For definiteness, hereafter we set J > 0; 
then the bound state exists for spin down. We denote the 
wave function of the bound state in the momentum rep- 
resentation as iy9k; the Schrodinger equation for it, which 
follows from the form of Hamiltonian mh, reads 



(eo - £k)'Pk = -JS^ 



<<5ki- 



(6) 



Clearly, the wave function must have the form (/Jk oc (eg ^ 
£k)~^- Substituting it into Eq. (^), we find the equation 
for the energy of the bound state (eo < 0), 



1 = JS 



1 



yielding 



So 



(27r)3 ek - eo 



-2C^A 



(7) 



(8) 



at ^ ^ 1. The power-law dependence, eo oc is a 
direct consequence of the excitations spectrum ek. The 
convergence of the integral in Eq. (^ at ^k ^ 
A validates the use of the simplified form of E'k. We 
emphasize that in order to obtain the bound state, it 
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was essential to go beyond the Born approximation.BLl 
The wave function of the bound state in the coordinate 
representation has the form 



that should be identified with 2\tij\. This way, we find 



sin kpr 



'€o/2C 



(9) 



where = vf/^ is the coherence length of a clean su- 
perconductor. 

The wave function of the localized state decays expo- 
nentially on the spatial scale which exceeds the coherence 
length by the factor Clearly, the effect of multiple 
magnetic impurities on the quasiparticle spectrum de- 
pends crucially on the characteristic distance between 
them compared to ^o/C- ■'^t lowest concentrations Us, 
separate magnetic impurities create a 5-peak in the quasi- 
particle DoS. With the increase of the concentration the 
peak broadens. As long as the average inter-impurity 
distance is large, the broadening is due to rare occur- 
rences of closely located impuritv-jpairs producing split 
levels due to their hybridization.o However, at higher 
concentrations, ?^s(Co/C)'^ ^ 1; the bound state on a 
given impurity overlaps with a large number of other lo- 
calized states and forms a well-defined impurity band. 
This large parameter allows one to treat the DoS in that 
band self-consistently. 



B. Impurity band 

We start with an estimate of the impurity band width 
W. A quasiparticle initially localized on a given impu- 
rity may hop on ~ ^t-sCCo/C)^ ^ 1 "nearest neighbors". 
Because of their large number, the return of the quasi- 
particle is improbable. Therefore, we may introduce the 
escape rate Ti for a given impurity. It will determine the 
level width, and we will associate it with the width of 
the band, ~ W. On the other hand, the escape rate 
from a given site is proportional to the DoS on the "re- 
ceiving" ones. The latter is inversely proportional to W. 
This way, the formula for the escape rate turns into a 
self-consistent equation for W (this is the essence of the 
self-consistent Born approximation for degenerate states, 
analyzed in detjyil .in the context of the broadening of high 
Landau levelsL3L3). To implement the scheme, we denote 
the hopping matrix element between two impurity states 
i and j, at positions Rj and R^, as t^ . The said self- 
consistency equation reads 



(10) 



(upon summing over positions j, the dependence on i dis- 
appears). The estimate of Uj should be obtained from a 
consideration of the two-impurity problem. Tunneling 
between two impurity sites splits the degenerate level 
in two. The resulting energy splitting depends on the 
inter- impurity distance Rij. Similarly to the conven- 
tional tight-binding problem, it is this energy splitting 



|(/j(R.y)|. To estimate the proportionality coef- 



ficient here, we note that the energy splitting reaches 



a value of the order of |£o| if kpRi 



1. Therefore, 



\tij\ ~ |£o| • l9'(R-y )l / |9'(1/^f)|- This estimate is con- 
firmed by a solution of the energy splitting problem, 
which is easily obtained by the generalization of Eqs. (^)- 
(0) to the case of two impurities.o 

Summation over the random positions of the impurities 



j results in ^ - \U 



2 

Remarkably, we used here only the value of the wave 
function ip{r) close to the impurity and the fact that 
is normalized. This is why one may expect that 
the obtained estimate of l^ijP remains valid in the 
presence of nonmagnetic disorder, as long as the corre- 
sponding mean free path is large compared to the elec- 
tron Fermi wave length. (This expectation is indeed con- 
firmed by the rigorous consideration of Sec . [Illl ) Sub- 
stituting the estimate of J2j l^ijP i'^to Eq. (|10D and re- 
placing kp/vp ^ t'o, we find expression for W up to a 
numerical factor. Together with this factor (which will 
be found later), the bandwidth in terms of measurable 
quantities Ug and Eq, reads 



21/4 



7ri/2 Wo A 



1/2 



1^0 1 



1/4 



(11) 



This expression is valid when the band is narrow, W <C 

ko|. 

The next question is the DoS i^b(£) for the Bogoliubov 
quasiparticles (hence the B subscript) inside the impurity 
band which can be found from the momentum-integrated 
Green function g{e) as 

VB{e) = Im.g(£ + iO), g{e) = / -— -3GB(k,£). 

(12) 

The Green function G_B(k, £) of the Bogoliubov quasipar- 
ticles is determined by the Hamiltonian, which now, in 
contrast to Hamiltonian (^), contains magnetic impuri- 
ties of finite concentration: 



k,(T ki ,k2 

(13) 

The Green function can be written in terms of the 
self-energy E(k, £) and we treat the self-epergy within 
the self-consistent T-matrix approximation,!^ which dia- 
grammatically amounts to the geometric series of Fig. |^: 



GB(k,£) 



S(k,£) = - 



1 



£ - £k - S(k,£) ' 

UgJS 



(14) 
(15) 



(the sign in front of J in the self-energy corresponds to 
quasiparticle spin down, which is of interest for us). In- 
tegrating Eq. ( p^ ) over momentum, we obtain the self- 
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FIG. 1. Self-energy E(k, e) corresponding to the self- 
consistent T-matrix approximation. The scattering on a sin- 
gle impurity (dashed lines) is treated in all orders of the per- 
turbation series. The self-consistency is achieved by inserting 
the bold lines for the Green functions, which means that they 
are calculated with the help of the presented self-energy. 



consistency condition for g\{s): 



JSgi 



n,JS 



(16) 



The impurity band corresponds to nonzero DoS, that is, 
to the nonzero imaginary part of g^, in some domain 
of negative e. Analyzing Eq. (|l^ in the limit rig <C 
Ci'oA (the Umit of narrow band), we find a semicircle 
DoS inside the impurity band: 



1^0 I 



W/2 



(17) 



which is centered at Eq, found earlier in Eq. (|^), and has 
the width given by Eq. (^ij) [note that it is the calculation 
leading to Eq. (^^ that establishes the numerical factor 
in Eq. (pT|)]. The total DoS inside the impurity band 
is equal to Ug, demonstratiag that each impurity brings 
one state into the system.tHi The same Eq. (|l6|) can be 
employed to demonstrate that the lower edge of the con- 
tinuous spectrum shifts from zero to Ec = ng/STr^j^QA at 
finite impurity concentrations, and the behavior of the 
DoS near the edge is 



(18) 



(The positive value of Ec is a manifestation of the level 
repulsion between the states of the continuum and the 
impurity states.) 

As we mentioned before, a well-defined impurity band 
exists [and can be described by the self-consistent T- 
matrix method leading to DoS (0)] if ns(Co/C)^ > 1, 
which translates into W ^ \eo\\/C{^/ Ep), where Ep \s 
the Fermi energy. With the increase of the impurity con- 
centration, the impurity band broadens, W oc y^s- 
Ug ~ C'^oA, the width of the band becomes of the same 
order as the distance from the center of the band to the 
continuum edge, resulting in the merger of the impurity 
band with the continuum of quasiparticle states. [Quan- 
titatively, the analysis of Eq. (|l^) shows that the energy 
gap between the impurity band and continuum exists if 
ns/(Ci/oA) < 167r/27 w 1.86.] The two conditions for the 
existence of the impurity band, ns{^Q/C,)'^ 3> 1 and iis < 
Ct'oA, are consistent since vq/^S^q/C,"^ ^ (Ep/AQ'^ >• 1. 
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FIG. 2. Different regimes for structure of the DoS due to 
polarized magnetic impurities of concentration ris . 



In the above discussion, we did not treat the self- 
consistency equation explicitly, assuming A that we used 
to be the self-consistent order parameter corresponding 
to the given impurity parameters. Actually, A is sup- 
pressed by magnetic impurities. A gapless state is re- 
alized at a critical concentration Ug ^ vqA/C,- At that 
concentration the average Zeeman energy associated with 
the magnetic impurities, 

E, = nsJS, (19) 

is comparable to A 00 Figure ^ illustrates the results 
of the above discussion. 

The above consideration of spin-polarized magnetic 
impurities can be easily generalized to the case of ran- 
domly oriented spins (see Appendix]^. It yields Eq. (A3) 
for the DoS — now summed over the quasiparticle spin 
directions — inside an impurity ban d of width W see 
Eq. (|a3)]. Comparing Eqs. ( |A3|) and ( |A4| ) with Eqs. (1^) 
and (|lip, wc notice the relation W — yV/-\/2 between the 
widths of the band. Indeed, randomness of impurity spin 
orientations leads to an additional cos(i?/2) factor in the 
hopping matrix element tij between two impurity states 
[see Eq. (|l^)] , and hence to the 1 /2 factor in the average 
{\tij\^), finally changing the width of the band by the 
l/\/2 factor. (Here is the angle between the spins of 
the two impurities.) 

Another difference of the random-spin case from the 
spin-polarized case is the absence of the average Zeeman 
energy. Upon increasing Ug, the gapless state is then 
reached due to the electron spin-flip scattering off mag- 
netic impurities. EJ The latter is characterized by a scat- 
tering rate 1/ts which, in the Born approximation, is 



1 

Tg 



2C'n, 

TTVo 



(20) 



The transition occurs when 1/ts = A, that is, at a con- 
centration of magnetic impurities ng ^ vq/S./Q'^. 



C. Effect of the impurity band on the 
electromagnetic field absorption 

At finite but low temperatures T, quasiparticles first 
populate the impurity band. This results in a finite dis- 
sipative part i7i(a;) of conductivity at arbitrary low fre- 
quencies^ We estimate cri(cij) at w ^ W^Jising a suitable 
versio nEa of the Mattis-BardeflB, theory.E3 We disregard 
possible spin selection rulesca'EJ by assuming that the 
spins of magnetic impurities are randomly oriented. 
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We may relate cti to the absorption power per unit 
volume, w = (l/2)(Tif^, where £ is the amplitude of the 
applied ac electric field. The absorption power can be 
evaluated as the product of the photon energy quantum 
UJ and the rate of creation of quasiparticlcs. The latter is 
oc £^ and is evaluated with the help of the Fermi golden 
rule. The perturbation leading to the excitation of the 
system is i? = (e/2)(Av + vA), where v is the operator 
of velocity of an electron and A is the vector potential of 
the ac field; its amplitude is A = £/oj. The absorption 
power associated with the photon-induced traiasitions of 
quasiparticlcs within the impurity band equalg23 



w 



2ttujB^ dEidE2[n{E2) - n{Ei)] 



X Vb{Ei)vb{E2){uiU2 + ViV2Y5{Ei 



Eo 



(21) 



Writing Eq. (|2l|), we accounted for the fact that the states 
in the impurity band do not correspond to a defined mo- 
mentum, and consequently the latter is not conserved in 
the transitions. To estimate the average over the states 
B^, we employ the arguments of the work of Mattis and 
Bardeen,c2l which allow one to express it in terms of an 
analog of the normal conductivity for the impurity band. 
That leads us to the estimate B"^ ^ £'^ (ev p j iS)^ (t j v) ^ 
where r is the relaxation time and i/ is a characteris- 
tic value of the DoS in the impuri ty b and. Substituting 
T -^XjW m^Av ^ Us/W [see Eq. (^], we find 



Tic V a; 



(22) 



In the case of "shallow" levels {C, <^ 1) we may also re- 
place the combination of the coherence factors, (miU2 + 
viV2)'^ ^ 1. At lowest frequencies, uj <^T, W, we may ex- 
pand the difference of the quasiparticle distribution func- 
tions in Eq. ( pl| ) as n{E2) — n{Ei) w —{dn/dE)uj. As a 
result, we find frequency-independent asymptotic forms 
for the dissipative conductivity at T <C and T 3> W, 
which we express in terms of the density riqp of quasipar- 
ticlcs in the impurity band [hence the (6) superscript]: 



^(6) _ f 4p^(ewf )V(rW^3)i/2^ at w < T < W^, 
'^^ ^ \ 4p^(ewF)V(rW^), at w < < T. 
In equilibrium, the density of quasiparticlcs is 



(23) 



'qp 



2 J dEvB{E)n{E) 



(24) 



{T/Wf^ g-(A-|eo|-W/2)/T^ y ^ 



rise 



-(A-|eol)/T 



at T » W". 



In the opposite limit of photon energies high compared 
to temperature, we may neglect n{Ei) in Eq. (^l|), finally 
finding 



exceeding the threshold 2A (each photon creates a pair 
of quasiparticlcs). The presence of states within the gap 
A reduces the threshold frequency and modifies the de- 
pendence of the dissipative part of conductivity cti on 
frequency for uj exceeding the threshold. Similarly to the 
above consideration, one may relate (Ti(cj) to the absorp- 
tion power (cf. Ref. p5|). 



w = 27rwS2 J dEidE2 

X iyB{Ei)iyBiE2){uiV2 - viU2fS{Ei + E2 - lj). (26) 



As follows from Eq. (A3), absorption is possible only if 
the frequency exceeds the threshold, 



W 

w > wth , wth = 2 ( A - |£o| - 



(27) 



To evaluate the absorption, however, the accepted ap- 
proximation of shallow levels is insufficient, because it 
yields (wif2 — ^1^2)^ — 0. As foll ows fr om the de- 
tailed calculation presented in Sec. Ill B 1 this factor 
{uiV2 — viU2)'^ ^ C^; the dissipative conductivity scales 
as a'"i\u!) (X {lu — Wth)^ near the threshold. 



III. DIRTY LIMIT 

We turn now to the practically important case of short 
elastic mean free time Tp caused by potential scattering of 
electrons and consider the dirty limit, A ^ I/t^. We will 
see that jjiatentiai disorder does not affect the electron 
spectrum.L^ However, its role is essential for transport 
properties (complex conductivity and impedance). 

In the dirty limit, magnetic impurities of small concen- 
tration caa. be described in the framework of the Usadel 
equationO Thaiorm of the Usadel equation, widely used 
in literature Mr ^1 corresponds to the Born limit for mag- 
netic scattering (with random positions and orientations 
of magnetic impurities). At the same time, the general- 
ization to arbitrary strength of magnetic scattering caa 
also be obtained, as shown by Marchetti and Simons.E§ 
Employing the standard 0-parametrization of the normal 
and anomalous quasiclassical Green functions. 



G = cos 9, F = sin^ 
we can write the generalized form asP^^^' 



sin 6* + A cos 9 - 



1 



sin 26* 



2ts 1 -I- C'' + 2C cos 29 



(28) 



= 0, 
(29) 



af\Lo) - n('')(ewi.)V(cjI^3)i/2^ T^io^W. (25) ^^^^^ ^ = ^pTp/^ is the diffusion constant and E is 



In an ideal superconductor at zero temperature, pho- 
ton absorption is possible only for photons with energy lu 



the electron energy measured from the Fermi level. We 
are interested in the homogeneous solution. Following 
Ref. 38, we define a "renormalized" energy and order 
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parameter, 



E 

A = A - 



cos 6* 



2ts 1 + + cos 261 ' 
1 sine* 
2^1 + C'' + 2C2 COS26'' 



(30) 
(31) 



With these quantities, the homogeneous Usadel equation 
becomes simple. 



ii; sin 6* + A cos 6* = 0, 
and can be formahy solved as 



cosf 



sm9 



(32) 



(33) 



where u = E/K. At the same time, this is not a solu- 
tion yet, but an alternative parametrization of the Green 
functions in terms of u. The Usadel equation actually 
determines the Green functions' dependence on energy 
E; hence, in Eq. ( p3| ) u itself depends on E. The lat- 
ter dependence is determined by the equation, directly 
following from Eqs. (^o|) and (^l|): 

M 1 + 7s 



E 
A' 



where 



1 



Is = 



(1+C2)2t,A 



£0 = 



i-e 



(34) 



(35) 



Here 7^ is the exact scattering rate due to magnetic im- 
purities, normalized by A, and eo can be shown to deter- 
mine the center of the impurity band, also in units of A. 
Since we focus on the case when each impurity is weak, 
C ^ 1, we simplify the above expressions as 



7s 



1 



£0 = 1- K 



(36) 



Once found, u can be immediately applied for analyz- 
ing the spectrum. The DoS per spin is given by 



^0 



= RcG{E) = Rc 



(37) 



It turns out that the potential scattering does not influ- 
ence th&^ergy spectrum, since Eqs. (^jk (p5|), and (|37)| 
coincideEa with those obtained by ShibaQ and RusinovQ 
in the clean limit (see also Appendix 



A. Green functions and the density of states 

Finite DoS corresponds to complex solutions of Eq. 
( p^ ) for u. Analyzing the behavior of the function in the. 
left-hand side (l.h.s.) of Eq. (H) at < u < 1, Shibafil 
found the domain of energies [the right-hand side (r.h.s.) 



of Eq. (^^] that cannot be matched by real values of 
u. At positive energy, this domain is the impurity band 
centered around Eq = |£o|A. Galculating u inside the 
narrow band, we find 



„IE) = ^ Re . 
for the DoS and 



1 - 



W = 4 



1 



21/4,^1/2 ^ly^A 



E — Eq 

W/2 



1/2 



at E>0, (38) 



£0 



1/4 



(39) 



for the width. The DoS at -E < is the mirror image of 
Eq. (^8|). The total number of states in the two impurity 
bands equals Ug . The DoS ( p8| ) coin cide s with the result 
derived in the clean limit [see Eq. (A3)] up to a factor 
of 4. This factor accounts for the difference in the def- 
initions of DoS. [The DoS (A3), defined for Bogoliubov 
quasiparticles existing only at _E > 0, corresponds to Eq. 
(|38| ) (electron component) folded together with its mirror 
image (hole component) and, additionally, summed over 
the spin projections.] 

As the concentration of magnetic impurities (and 
hence 7^) grows, the upper edge of the band merges with 
the continuum at ns/(Ci^oA) = 167r(2Av/3- 1)^/2 « 3.O6 
[this result can be obtained from Eq. (|4|) and differs only 
numerically from the condition obtained for the case of 
polarized magnetic impurities in Sec. [IB|. At higher con- 
centrations the AG regime is realized, with the gap width 
reduced due to the scattering off magnetic impurities. As 
follows from Eq. (|3^ ) , the quasiparticle continuum starts 
at energy lower than the single-impurity bound-state en- 
ergy Eq, and the gap Eg is reduced progressively with 



the increase of 7^ . 
with — 1 in Eq 



At7s/'»c^ 



one may replace — Cg 



m. Then it simplifies to the form 



l-7s 



1 



E 
A' 



(40) 



exactly as considered by Abrikosov and Gor'kov.0 Note 
that 7s itself can still be small, and we mainly focus on 
that case: 



7s < 1- 



(41) 



The AG regime, corresponding to a single gap in the 
spectrum, has two characteristic features, essential for 
the low- frequency dissipation: first, the gap is suppressed 
(the relative scale is determined by 7s ^^), and second, 
the square-root BCS singularity of the. DoS at = A is 
smeared. The gap Eg in units of A isEl 



""^-(1-7?/^ '-l-k^^ 



A V - J - 2- • 
Solving the Usadel equation at energiesclose to the ac- 
tual gap, {E - Eg) < 7s^^A, we obtairEl 



GiE) 



1 l2{E^Eg) 



1/3 



2/3 



3A 



(43) 
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The square-root behavior of the DoS [see Eq. (pTI) ] in the 
vicinity of the gap was found in Refs. ^ and |45| . 

At low concentration, magnetic impurities substan- 
tially modify the spectrum in a relatively narrow domain 
of energies, \E — A\ < 7s^"^A. We take i? = A as a 
reference point within the domain, and find G and its 
derivative with respect to energy: 



G(A) 



1 



G'(A) = 



37s 3 • 24/37, 



1/3- 



(44) 
(45) 



The maximum of the DoS is achieved slightly above 
E ^ A. At the same time, considering the derivative 
G'(A) and the BCS solution that makes it possible to 
approach the maximum from the side of higher energies, 
we conclude that the maximal DoS is of the same order 
as the value determined by Eq. (H). This result was 



found in Ref. 25 



The results for the DoS are summarized in Fig. Note 
that the square-root behavior of the DoS near the gap 
edge [following from Eq. (^)], being applied at ii^ = A 

(which is already beyond its applicability range), yields 
— 1/3 

i'{A) = , which only slightly differs from the accu- 

rate result i^A) = (^3/24/3)77^/^ « 0.77^^^^ following 



from Eq. (|44|). 

In the evaluation of dissipation with the help of Eq. 
(p7|), a good approximation for the DoS is given by 



1 

2/3 



3A 



Re G{E) 



V3 



E 



Eg<E<A, 
A<S< A(l+7,'/'), 
A{l + -f^s^^)<E. 



(46) 

At the matching points, the above expressions match by 
the order of magnitude, being different only by numerical 
factors close to 1. Therefore, this is a rather accurate 
approximation for calculating the dissipation up to order- 
of-one numerical factors (at the same time, in the limiting 
cases that we consider below, the numerical factors will 
be asymptotically exact). 

At small and 7s, magnetic impurities significantly 
affect the electron spectrum only at energies close to A. 
In that energy range [corresponding to Eq. ( p^ ) and to 
the first two lines in Eq. (^) ], w e have u ~ I and hence 
F K iG [as follows from Eq. (p3|)], so that finally ImF sa 
ReG. Outside that domain, the functions F and G in 
the leading order have the conventional BCS form [with 
ImF(£;) « A//B2— ^ instead of the last line in Eq. 
(!)]■ 

Note that we have discussed ReG(£') and ImF(£') at 
-E > 0. At the same time. Re G is even while Im F is odd 
with respect to E. 




i-(3/2)r 



FIG. 3. Schematic (not numerically exact) DoS for a super- 
conductor in the AG regime with a small concentration of 
magnetic impurities. The bold sectors of the curve are known 
analytically. 



B. Dissipative conductivity 

Having found the Green functions, we are ready to cal- 
culate the dissipation. The expression for the response 
kernel of the current to the external electromagnetic field 
is derived in Appendix From the result for the kernel, 
we extract the expression for the dissipative part of con- 
ductivity: 



- I dE\ tanh 



2T 



tanh ■ 



E_ 

2r 



X [ReG(£;+)ReG(£;_) + ImF(£;+)ImF(£;_)] , (47) 

where E± — E ± ijj/2 and = 2e^VQD is the Drude 
conductivity. This formula is the dirty-limit counterpart 
of Eqs. (^ ) and (26) for the absorption power in the clean 
limit. The combination of the Green functions G and F 
here replaces the product of the quasiparticle densities of 
states and the combinations of the coherence factors in 
Eqs. (|l|) and (|2|). 

The analysis of the present section is be based on Eq. 
(^). At the same time, the results can be easily con- 
verted into the l angua ge of the surface impedance, as 
discussed in Sec. HID. Small but finite dissipative con- 



ductivity implies a small but finite real part of the surface 
impedance. 



1. T — 0: Transitions across the gap 

At zero temperature the dissipation is absent below a 
certain threshold frequency. In the absence of magnetic 
impurities, uj must exceed 2 A so that a Cooper pair can 
be broken and two quasiparticles are created above the 
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gap. At frequencies slightly above the threshaLd_(tJ 
2 A) ^ A, the dissipative part of conductivity isHElcJ 



2 V2A 



(48) 



In the presence of magnetic impurities, the BCS sin- 
gularity of the Green functions is smeared out and hence 
the threshold behavior of the dissipation changes. The 
threshold frequency in the AG regime is 2E q, a nd the 
nonzero contribution to the integral in Eq. ( |47|) comes 
only from the region of energies < E < {uj/2 — Eg), 
where the arguments of the Green functions are close to 
iEg SO that, at first sight, we can employ the square- root 
form for ReG and ImF [the first line in Eq. (^)] and 
take into account that ReG(i<^) is even while Imi^(_E) 
is odd. However, approximating ReG(£^) — lmF{E) at 
-E > we find that the integrand in Eq. (^7|) vanishes in 
the main order, so we should take into account a small 
difference between ReG and Imi^. This is easily done 
on the basis of the normalization condition + — 1. 
Note that it is the large absolute value of G near Eg [see 
Eq. (H)] that leads to F w iG and hence to Im F sa Re G. 
At the same time, the normalization condition immedi- 
ately gives also a correction to this result: 



ImF « ReG 



^^2G- 



(49) 



Then we find nonzero dissipation which is quadratic in 
deviation from the threshold at (a; — 2Eg) <^ 7s^'^A: 



67; 



2/3 V 2E, 



(50) 



This result was obtain in Ref. ^ The quadratic thresh- 
old behavior of ai is a consequence of the square-root 
threshold behavior of the DoS, and experimental results 
of Ref. H seem to confirm this (Ti(aj) dependence. 

At lowest concentrations of magnetic impurities, the 
impurity band exists inside the gap, so that the thresh- 
old frequency for absorption is determined by transitions 
in which a Cooper pair is broken and two quasiparticles 
are created inside the impurity band. The threshold fre- 
quency for this process is Wth — '^{Eg — W/2), which is Eq. 
( 2^ written in different notation. Just above the thresh- 
old, at {uj — Wth) ^ W, only the edge behavior (near an 
edge of the impurity band) of Re G and Im F enters Eq. 
(p7|). According to Eq. (^), both functions near an edge 
grow as a square root. We again find a cancellation in the 
main order and have to employ Eq. (^). The functional 
form is thus the same as in the AG regime, and the result 
differs from Eq. (E^) only by coefficients: 



2Cn 



2^2 



Co 



(w - Wth) 



(51) 



Next we consider the case of finite (but low) tempera- 
tures, when the threshold for absorption is absent since 
equilibrium thermally excited quasiparticles can absorb 
photons of arbitrarily small energy. 



2. Low T: Transitions within the AG states 

We start with the consideration of the AG regime in 
the limit of small 7s, so that quasiparticle states exist 
only above the gap Eg, which is weakly suppressed com- 
pared to A, as described by Eqs. ( ^ ) and (^). Most 
of the results of this subsection were obtained in previ- 
ous publications, although sometimes in the framework 
of physically different problems and with the help of dif- 
ferent techniques. Our consideration makes it possible to 
incorporate them in a unified manner. 

In the limit that we consider, ReG(£') and ImF(i?) 
entering Eq. (|4^) are nonzero only above Eg which is 
close to A, so that at low frequencies and temperatures. 



a;,r< A, 



(52) 



two simplifications are possible: (a) we can substitute 
the tanh's with their asymptotic exponential forms, and 
(b) since the exponentials limit the integration to the 
region of |F — A| < T, where lmF{E) is approximately 
equal to ReG(i?), the contribution of the (ImFImF) 
term is the same as the contribution from the (Re G Re G) 
term. Therefore, we can simplify Eq. (|47|) as 



CTO 



= ^jdE (e-^'^ - e-(^+")/^) 

y.ReG{E)KeG{E + uj). (53) 

At the same time, at cj ^ A the external electromagnetic 
field cannot excite quasiparticles across the gap; there- 
fore, the dissipation will be due to quasiparticles that are 
already excited due to low but finite temperature. It is 
thus natural to express the results for the dissipation in 
terms of the equilibrium quasiparticle density, 
00 

nqp = 41/0 j dEReG{E)-^j^-^, (54) 


for which under the same assumptions (T <C A, 7^ <C 1) 
we obtain two limiting cases: 




'"qp _ 

I'D 

1^0 



at T < 72/3A, (55) 
at T > 7s/^A. (56) 



Below we consider various relations between the three 
energy scales uj, T, and 7s ^'^ A (while all of them are 
much smaller than A). Starting with the case of rela- 

2/3 

tively strong magnetic scattering, T ^ 7s A, we note 
that, due to the exponentials, the main contribution to 
the integral in Eq. ( p3[ ) comes from the narrow region 
above Eg, where the square-root asymptote for the DoS 
is valid [see the first line in Eq. (^)]. As a result, 

ai{uj) 8T 



.-^./T^inhf^Ui^-^ 



\2TJ 



\2TI 



(57) 
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where Ki is the Macdonald function (the modified Bessel 
function of the second kind). Asymptotic forms of the 
Macdonald function at large and small values of the ar- 
gument describe the two limiting cases for the ui/T ratio, 

so that finallyi-a, 
result of Maki,yJi 



2/3 

^ T ^ 7s A we reproduce the 



0-1 (g;) 




while at T < w < 7s A we find 



2/3 

7/ 



A 



1/2 



J^oA' 



(58) 



I/O A' 



(59) 



The result (57) was actually obtained by Larkin and 
Ovchinnikov iu-A physically different but formally equiv- 
alent problem.E^ They considered a superconductor with 
inhomogeneous effective interaction between electrons. 
For small-scale inhomogeneities the problem turns out to 
be mathematically equivalent to the case of homogeneous 
effective interaction in the presence of magnetic impuri- 
ties. The results of Ref. ^ thus apply to the problem 
considered in the present paper. 

In the limit of weak magnetic scattering, uj,T 3> 
7s^^A, the main contribution to the integral in Eq. ( |5^ ) 
comes from the region of energies above the smeared BCS 
singularity; in that region we can apply the unperturbed 
result for the DoS [see the last line in Eq. (p6|)]. We thus 
reproduce. the res ults obtained in the absence of magnetic 
impuritiescllE3ca (we have^o take into account the corre- 
spondence noted by NamElbetween the extreme anoma- 
lous limit in the clean caseEd and the dirty case to which 
our calculations refer): 



(60) 



Asymptotic forms of the Macdonald function Kq at large 
and small values of the argument describe the two limit- 
ing cases for the u/T ratio, so that finally at 7s^'^A <C 
w <C T we find 

Tijoj) 1 /A 
(TO \T 



3/2 



1^0 A '-fu! 



(61) 



(where 7 « 1.78 is Euler's constant), while at 7, 
T <^u) the result is 



2/3 



3/2 



I/O A 



A « 



(62) 



The logarithmic divergence in Eq. ( |6l] ) at cj — > is due 
to the overlap of two square-root BCS singularities of the 
DoS. 

Similarly to the case of Eq. (E32 
Lo the second exponential in Eq. 



can be neglected 
while the w-shifted DoS is nearly constant, TieG{E + 



at T < 7s/^A < 



result (|62|) in terms of n. 
itself are different 



at 



js^^A < w [Eqs. (H) and 



Lo) « a/ A/2a;, in the essential region of integration. As 
a result, the remaining integral in both cases is exactly 
the same as the one determining the quasiparticle density 
(^J) at T ^ Eg. Therefore, we obtain exactly the same 
qp, although the expressions for 

-fs^^A < T < w and T < 

]55| ), respectively]. This result 
actually follows from Ref. SqT 

Finally, at oj ^ 7s^'^A <^ T, relatively large temper- 
ature allows rather wide region of energy integration in 
Eq. ([5^), while the shift of the two DoS can be neglected. 
Similarly to the case of Eq. (|6^) , a logarithmic singular- 
ity appears in the integral; however, now it is cut off not 
due to the finite shift w [as it was the case in Eq. (|6ll)] 
but due to smearing of the square-root BCS singularity 
of the DoS by magnetic impurities. As a result, with the 
logarithmic accuracy we find 



o"! (g;) 
c^o 



1 



/2ti 



3/2 



"qp 
I/O A 



In 



T 



7s A 



(63) 



The results for aiiuj) at cj ^ A are valid as long as the 
quasiparticle population may be described by the Boltz- 
mann distribution with an effective temperature. The 
quasiparticle chemical potential should not necessarily be 
and the temperature T may deviate from equilibrium 
as long as T <C A. 



3. Low T: Transitions within the impurity band 

At equilibriu m, the results for the dissipation ai cal- 
culated in Sec. Ill B 2 for the AG regime in the limit of 
low temperatures always contain an exponentially small 
factor exp{—Eg/T), since the absorption of the electro- 
magnetic field is due to the quasiparticles thermally ex- 
cited above the gap. At the same time, as we discussed 
above, at low concentrations of magnetic impurities a 
band of quasiparticles states appears below the contin- 
uum. This means that at temperatures much smaller 
than the distance between the band and the continuum, 
the dominant contribution to the dissipation will be due 
to quasiparticles residing inside the impurity band. Be- 
low we consider this case, which is the main focus of our 
interest. 

Our aim now is to calculate the dissipation (^^ for 
the case when it is due to transitions inside the impurity 
band; hence, uj < W. At the same time, we assume the 
band to be narrow, so that its width, which can be writ- 
ten as W — 4-^C7sA, is much smaller than the distance 
from the center of the band Eq to the bottom of the con- 
tinuum A: Since Eq — {1 — 2C^)A, the condition can be 
written as 



il^' « c 



3/2 



(64) 



The contribution of the impurity band is dominant at 
T C'^A, when the continuum contribution is exponen- 
tially smaller due to the exp(— 2C^A/r) factor. 
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Since we consider energies of the order of Eq, that is, 
of the order of A, and ReG{E) — Im F {E) inside the 
narrow band (as discussed in Sec. HI A), we can make 
the same simpHfications that lead us to Eq. (|5^). The 
only difference is that now the region of the energy inte- 
gration is confined to within the impurity band. Taking 
into account the explicit form of the DoS inside the band 
[Eq. (^8|)] and making the variable of integration dimen- 
sionless, we obtain 



while at T ^ Lo the result is 



(b) 



(71) 



At {W — uj) <C r (while W and uj are not necessarily 
close to each other), the exponential in Eq. (^5|) can be 
replaced by unity in the region of integration, and the 
remaining integral can be written as 



















^ / 


de e~'^^\j 


-(1-w) 





_iEo . , w 
e -r smh — 



1 - 



UJ 



(65) 



The integral contains the product of the two semicircle 
DoS (|3^) shifted by the frequency uj. The overlap exists 
and produces nonz ero diss ipation only aX uj <W . 

Similarly to Sec. |III B % we express the final results in 
terms of the equilibrium quas iparticle density (within the 
band). Calculating Eq. (54) for the impurity band, we 
find 



(b) _ -Bo/Tr 



2T j 



(66) 



where Ii is the modified Bessel function of the first kind. 
In the two limiting cases for the W/T ratio, we obtain 



"qp 



3/2 



Eq-W/2 



at T < VF, (67) 
at T > M^. (68) 



Below we consider various relations between the three 
energy scales, w, T, and W (while all of them are much 
smaller than A). Starting with the case of relatively large 
band width, a;,T ^ W, we note that, due to the expo- 
nentials, the main contribution to the integral in Eq. ( |65| ) 
comes from the narrow region above the lower limit of in- 
tegration where the simplified square-root result for the 
DoS is vahd: v{E) cx ^/E~ {Eq - WJ2). This situation 
is very similar to the case of Eq. (B^ (cf. the discus- 
sion above that formula). As a result, we obtain the 
same functional form [with a different prefactor and with 
{Eq — W/2) instead of Eg in the exponent]: 



(w) Un^T 



-\V/2 



r2„2 



M/3 



(69) 

Asymptotic forms of the Macdonald function at large and 
small values of the argument describe the two limiting 
cases for the uj/T ratio, so that finally at a; ^ T <C 
we find 



(c) 



(b) 



7r3/2 j,2^3/2yl/2 



(70) 



(cj) 16n^^{W + uj) 



Stt^i^^W^uj 



sinh — 
2T 



(72) 



Here K and E are the complete elliptic integrals of the 
first and second kind, respectively, and their argument is 



k = 



W - 



W ■ 



(73) 



In the limit of small frequency and high temperature, 
a; < < T, with the help of Eq. (||), expression (|7|) 
can be simplified as 



{uj) 



Ab) 



Co 



37r2 vlWT' 



(74) 



There is an upper threshold for impurity band absorp- 
tion at a; = W . Close to that threshold, Eq. ( [7^ ) gives 

a^^\uj) cx iW — ujY9{W — uj), where is the Heaviside 
step function. 

Comparing Eqs. (7^, (|7l|), and ( [t^ in the dirty case 
with Eqs. (p3) and (25h in the clean case, we notice that 
the dependences of the dissipative part of conductivity 
on UJ, T, and W are the same. At the same time, the 
overall coefficients in the two cases are difi^erent due to 
different limits with respect to potential scattering. 



C. Magnetic impurities as traps for nonequilibrium 
quasiparticles 

At low temperatures, the density of equilibrium quasi- 
particles becomes negligible, and in reality riqp isrdpm- 
inated by extraneously produced quasiparticles .e3 In 
quantum devices based on Josephson junctions, the ad- 
verse effect of quasiparticles on the coherence may be mit- 
igated by the inclusion of traps — sections of the super- 
conducting leads with reduced value of the gap. We may 
pose a similar question with respect to the impedance: 
Can magnetic impurities reduce the dissipation caused by 
nonequilibrium quasiparticles? Such a possibility is most 
intriguing in the range of frequencies w < |£o| [seeEq. (H)] 
and quasiparticle densities nqp <^ C'^oA (see Fig. ||). In 
the following, we use the results of the two previous sub- 
sections to elucidate the dependence of the normalized 
dissipative conductivity ai/ao on n^. 
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For estimates, we assume the density of quasiparti- 
cles riqp fixed, and their effective temperature low, T <^ 
|£o|/ lii(C^'o^/'i'qp)i so that thermal ionization to the con- 
tinuum does not prevent magnetic impurities from trap- 
ping the quasiparticles. In the absence of magnetic im- 
purities, the dissipative response is defined by Eq. (62). 
At smah density, impurities reduce by the density of 
quasiparticles at the bottom of the continuous spectrum, 
yielding cti/ctq ~ (l/y2)(A/a;)3/2(nqp - n,)/(iyoA). In 
the interval of densities < ^ C'^oA, the impu- 
rity band may accommodate all the quasiparticles, while 
being too narrow to allow absorption of the field of fre- 
quency uj. This low- absorption regime stretches up to 
the densities Ug < (7r/8)(Ci'oA)(ci;/|£o|)^. At higher Us, 
the impurity band is broad enough to allow for intra- 
band absorption; compared to the benchmark Us = 
dissipation level, Eq. ([Tl]) yields (Ti/cfq smaller by a fac- 
tor ~ (w/|eo|)("-s/C'^oA)^/^. Upon further increase of TT-s , 
the impurity band merges with the continuum. Deep in 
the AG regime, Us ^ C'^oA, one may use Eq. ( p^ for 
the estimate of dissipation; compared to the Ug = dis- 
sipation level, Eq. ( p9[ ) yields cri/uo smaller by a factor 
~ (w/|eo|)(Ci^oA/n,f73_ 

To summarize the above estimates, in the selected re- 
gion of frequencies, cu < |eo|i and nonequilibrium quasi- 
particle densities, Tiqp <^ C^oA, magnetic impurities in- 
deed reduce the dissipative part of conductivity due to 
the reduction of the DoS available for the quasiparticle 
transitions. 



D. Surface impedance 

Finally, we reformulate the above results for the 
complex conductivity cr(w) in terms of another impor- 
tant experimentally measurable quantity, the surface 
impedance:OEj 



Z(^) 



/ Attuj 



(75) 



The impedance can be written as Z = i? — iX , where R 
and X are the surface resistance and reactance, respec- 
tively. The complex conductivity a = ai + i<T2 is given 
by Eqs. ( pT6| ) and ( pT^ ). 

Since we always consider the case of weak dissipation, 
ci <J2, we also obtain R X. The response of the 
superconductor is then almost purely reactive with 



Xiio) = 



(72 (w)c2 



(76) 



Under assumptions of Sec. Ill B , we can calculate ct2{(jj) 
in the simplest zero-temperature BCS model without 
magnetic scattering (since corrections due to low tem- 
perature and weak magnetic scattering, T, 7s^"^A <^ A, 
are negligible). At the same time, the new qualitative 



feature arising due to finite ai (calculated in Sec. Ill] 
is the dissipation, i.e., small but nonzero R: 



XHu;y 



(77) 



In the near-threshold situation of Sec_ [II B 1 , we are 
interested in frequencies w « 2A, thenEHl 172 (2 A) ~ ctq, 
and we find 



, , V27rA , , 



ca. 



(78) 



IIIB2 



I n the sm all-frequency limit pConsidered in Sees 
and |III B ^ , at i:j < A, we findE3 0-2 (w) = (7rA/a;)CTo and 
finally obtain 



R{lo) 



:(aoA)3/2 



(71 (w). 



(79) 



Thus, with the help of Eqs. (|7|) and (^), the results 
for (7i ( uj) (th e dissipative part of conductivity) calculated 
in Sec. IIIB are directly translated into the surface resis- 
tance. 



IV. CONCLUSIONS 

We have employed the quasiclassical approach to study 
dissipation in superconductors due to small concentra- 
tion of magnetic impurities. The superconductor was 
assumed to be in the dirty limit with respect to the po- 
tential (nonmagnetic) scattering. The concentration of 
magnetic impurities was assumed to be small enough so 
that the gap suppression is weak. Empkaying the exten- 
sion (proposed by Marchetti and SimonsCS) of the Usadel 
equation beyond the Born approximation for magnetic 
scattering, we have considered in a unified manner both 
the AG regime (with a continuum of states above a gap) 
and the limit of lowest magnetic impurities' concentra- 
tions where the impurity band due to overlap of localized 
impurity states is formed below the edge of the contin- 
uum. The deviation from the Born limit was assumed to 
be finite but small, so that the subgap impurity states lie 
close to the continuum edge. 

Our main results refer to the limit where the tempera- 
ture T, the frequency uj, and the gap suppression 7s^'^A 
are all much smaller than the BCS gap A in the absence 

of magnetic impurities. At the same time, the relation 

2/3 

between T, uj, and 7s A was assumed to be arbitrary, 
and we have obtained explicit analytical expressions in 
various limiting cases. 

Our results can be expressed in terms of quantities de- 
scribing the response of superconductor to external elec- 
tromagnetic field: the dissipative conductivity cti (the 
real part of the complex conductivity a) and the surface 
resistance R (the real part of the surface impedance Z). 

In the limit of small temperatures and small frequen- 
cies that we considered, at equilibrium, the dissipation is 



12 



always proportional to the density of thermally excited 
quasiparticles and thus exponentially suppressed due to 
gapped character of the spectrum. On the other hand, 
fluctuations in positions of magnetic im purit ies can lead 
to finite DoS below the mean-field gap£30 Investigat- 
ing the effect of the "tail" states on the dissipation is 
beyond the scope of our work. We only note that since 
the number of tail states is small (compared to the num- 
ber of the "mean-field" states considered in this paper), 
there is a wide region of applicability for our results at 
subgap temperatures T/A <C 1, excluding only ultra-low 
temperatures. 

Note added. During preparation of this paper, we be- 
came aware of a preprint by Kharitonov et alE3 address- 
ing the same topic with a complementary approach. In 
that work the surface impedance of superconductors with 
magnetic impurities is evaluated numerically, with em- 
phasis on the superconducting gapless regime having low 
energy states in the electron spectrum. 
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Appendix A: Nonmagnetic disorder does not affect 
the electron spectrum 



1. Magnetic impurities with randomly oriented 
spins 

In the clean limit, the Hamiltonian describing quasi- 
particles near the edge of the spectrum and interacting 
with randomly oriented weak magnetic impurities is 



+ JS J2 «L(wnJ■)akV'e^^'''-'')^^ (Al) 

k,k' ,a",cr' .j 

where s is the vector of the Pauli matrices acting in the 
spin space and is the (random) orientation of the spin 
of impurity j. [The Hamiltonian (|l^) for polarized mag- 
netic impurities corresponds to all = z.] The T-matrix 
approximation for this model yields the self-consistency 
equation 



JSg = - 



\ 



|£o| 



(A2) 



for the momentum-integrated Green function .g(e), which 
is now diagonal in spin space. At low concentration of 
weak magnetic impurities (C ^ 1 and 7^ ^ 1) and for 
energies close to the ed ge o f the continuous spectrum 
(i.e., |u — 1| ^ 1), Eq. ( [A^ ) agrees with Eq. ( p^ that 
was derived in the dirty limit provided one identifies 
g ^ —2i'Kh'Q cos 9 [se e also Eq. (^3|)]. In particular, the 
analysis of Eq. (A2) in the limit Ug ^ C^'oA yields a 
narrow impurity band with 



£0 



W/2 



(A3) 



for the total DoS (summed over quasiparticle spins) and 
1 



W = 4- 



21/4^1/2 \^j/oAy vi£o 
for the width [which coincides with Eq. (M 



£0 



(A4) 



In this Appendix we demonstrate that the electron 
spectrum in the presence of magnetic impurities, regard- 
less of their polarization, is insensitive to the nonmag- 
netic disorder. To see thi s we consider separately two 
different models. In Sec. A 1 we derive explicitly the 



equation determining the DoS in clean superconductors 
with randomly oriented spins and we establish the di- 
rect relation with the result obtained in Sec. [II (dirty 
superconductors, randomly oriented spins). In Sec. |A 2| 
we formulate the Usadel equation (dirty limit) for super- 
conductors with polarized magnetic impurities and find 
that the DoS corresponds to the one obtained in Sec. || 
(clean limit, polarized magnetic impurities). 



2. Spin-polarized magnetic impurities 

Following the derivation of Eq. ( p9| ) by Marchetti and 
Simons,K3 we find that the Usadel equation for the case 
of polarized magnetic impurities is 



iaEz sin 0„ 



1 - + SictC cos 9a 



iEsm9„ + Acos9^-z Z — — = 0. 

(A5) 

Here 9^ parametrizes the quasiclassical Green functions 
that are different for electrons with spins up and spins 
down. At low concentration of weak magnetic impuri- 
ties (C <C 1 and 7s ^ 1) and for energies close to the 
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edge of the continuous spectrum, Eq. ( |A5| ) in a uniform 
superconducting state can be rewritten as 



JSg„ 



._nsJS_ 



(A6) 



in terms of the spin-resolved momentum-integrated 
Green functions = —2i'KVQCOs9a- The equation for 
down spin is the same as Eq. ( [l6| ) derived in the clean 
limit, while the one for spin up does not produce an im- 
purity band. 



Appendix B: Derivation of the response kernel 

The Usadel equation in a superconductor reads 

Dd(gdg) - hdtg-dt'gfs - [A,. 9] - = 0, (Bl) 





A = 



A* 



(B2) 



where the quasiclassical Green's function g{r,t,t') is a 
matrix in the Nambu (Pauli matrices Ti) and Keldysh 
spaces (Pauli matrices &{), 



gK 

g^ 



(B3) 



Keldysh 



and the self-energy describes magnetic scattering.E^I 
The covariant derivative c? = V — i(e/c)[f3A, .] depends 
on the vector potential A(r, t) (we choose the gauge with- 
out the scalar potential) . Time convolution is implicit in 
the matrix product. The superconducting gap A(r, t) en- 
ters the matrix (in the Nambu space) A and solves the 
self-consistency equation 

A{r,t)^-^^tr{a,f^g{r,t,t)), (B4) 
where A is the BCS pairing constant. 



We look for the solution to Eq. (Bl) perturbatively in 



A, expanding g = +g^^^ + ... and A = A'") -I- A^^^ 
... In the zeroth order, g^^^ describes an equilibrium uni- 
form superconductor. 
In the next order. 



a.g = V.9W-^[f3A,.g(«)], 



(B5) 



so that 



d{gdg) = 9^'^ V2g(i) - ^ [f3 VA, j (B6) 

in the same order. If we choose the London gauge with 
VA = (accompanied by the requirement of nA = at 
the surface), then the Usadel equation in the first order 
reads 



£gW = [AW,5(")] 



(B7) 



The l.h.s. of this equation is a linear operator £ acting on 
g^-^\ while the source term in the r.h.s. is proportional to 
A^'^K Taking the first order also in the self-consistency 
equation (B4), we see that our choice of the gauge leads 
to A^^) = and g'^^ = 0. (It is-actually well known that 
A^^) = in the London gaugeE3) 

The current flowing in the superconductor reads 

j(r,i) = ^tr(aif35C?5)r,,^, (B8) 

(where ctq = 2e^VQD is the Drude conductivity), which 
in the lowest order in A reduces to 



(0) 



j(r,i) = -^tr (aif35(°)T3Ag 



Making the Fourier transform, we obtain 
j(k,c.) = -ig(c.)A(k,w), 



(B9) 



(BIO) 



with the response kernel [we omit the (0) superscript of 
the Green function for brevity] 



Q{uj) ^'-^JdEtr [a,T,g{E)r39{E -u)]^'-^ j dEtr [f^g" {E)T^g'' {E - uj) + f^g'' {E)f^g\E - 00)] . (BU) 

Note that independence of Q on k implies locality of relation between j and A in the dirty limit. 

The Keldysh component of the Green function can be expressed in terms of the retarded and advanced ones, while 
the advanced component can be written in terms of the retarded one, 

g'^iE) = {g"{E) - g^E)) tanh g^ = -f35^^f3, (B12) 

so that finally only the retarded component enters. In the case without the superconducting phase, it has the form 



5« = Gf3 -f Ffi 



(B13) 
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[G and F are the functions introduced in Sec. [II, see Eq. (p8[)] and we obtain 

G{E+)RcG{E-) - iF{E+)lmF{E-] 



Q(uj) i f , ( , 

= - dEi tanh 

ao 2 7 1 2r 



tanh ■ 



2T 



G*{E-)ReG{E+) +iF*{E-)linF{E+] 



where E± — E ± uj/2 [note that the integrand in Eq. 
(B14) is even with respect to E]. The difference of this 
expression for the kernel from the resuhs of Refs. p0|-|2^ 
is that our expression is written in terms of the quasi- 
classical Green functions (in the dirty limit) that can be 
directly found from the Usadel equation. A similar ex- 
pression was employed in Ref. although it was for- 
mulated in terms of contour integrals in the plane of 
complex E and the cuts in the plane were made from 
±Eg, the values of the AG gap; thus, the AG regime 
with a single gap was explicitly assumed. Our expression 
(B14) is more general, applicable also in the case when 
the impurity band is present and in the case of gapless 



(B14) 



r 



superconductivity. 

Defining the complex conductivityEHI a = ai + ia2 that 
determines the (local) response of current to electric field. 



j(k,L^) =a(L^)E(k,w) 



we see that 



iQ{uj) 



(B15) 



(B16) 



Dissipation is determined by the real part of conductivity, 
ai{u] ) — R eg(aj) = — Im(5(w)/tj, for which, according to 
Eq. pll) , we obtain Eq. (g^) given in the beginning of 
SecPIBl 
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